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Abstract. Using modern nucleon-nucleon interactions in the description of 
the A = 3, 4 nuclear systems the x 2 per datum results to be much bigger than 
one. In particular it is not possible to reproduce the three- and four-nucleon 
binding energies and the n — d scattering length simultaneously. This is one 
manifestation of the necessity of including a three-nucleon force in the nuclear 
Hamiltonian. In this paper we perform an analysis of some, widely used, three- 
nucleon force models. We analyze their capability to describe the aforemen- 
tioned quantities and, to improve their description, we propose modifications 
in the parametrization of the models. The effects of these new parametrization 
are studied in some polarization observables at low energies. 



1 Introduction 

Realistic nucleon-nucleon (NN) potentials reproduce the experimental NN scat- 
tering data up to energies of 350 MeV with a x 2 per datum close to one. However, 
the use of these potentials in the description of the three- and four-nucleon bound 
and scattering states gives a x 2 P er datum much bigger than one (see for example 
Ref. p]). In order to reproduce correctly the three-nucleon bound state energy, 
different three-nucleon force (TNF) models have been introduced as the Tucson- 
Melbourne (TM), Brazil (BR) and the Urbana IX (URIX) models [21 Si]. These 
models are based on the exchange mechanism of two pions between three nucle- 
ons. In the case of the TM model, it has been revisited within a chiral symmetry 
approach [5], and it has been demonstrated that the contact term present in 
it should be dropped. This new TM potential, known as TM', has been subse- 
quently readjusted [6] and the final operatorial structure coincides with that one 
given in the TNF of Brazil. TNF models based on irp and pp meson exchange 
mechanisms have also been derived [7] and their effects have been studied in 
the triton binding energy [Hj. More recently, TNFs have been derived [9] using 
a chiral effective field theory at next-to-next-to- leading order. A local version of 
these interactions (hereafter referred as N2LOL) can be found in Ref. |10j . At 
next-to-next-to-leading order, the TNF has two unknown constants that have 
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to be determined. It is a common practice to determine these parameters from 
the three- and four-nucleon binding energies (i?( 3 H) and -B( 4 He), respectively). 
It should be noticed that in this procedure, the three- and four- nucleon sys- 
tems are described in the framework of the non relativistic quantum mechanics. 
Relativistic corrections to the few-nucleon binding energies have been studied in 
Ref. and, recently, the three-nucleon Faddeev equations have been solved in 
a Poincare invariant model [12]. These efforts are directed to establish if some of 
the discrepancies observed between experimental data and theoretical descrip- 
tions, as for example the minimum of the N — d differential cross section, can be 
reduced if relativistic corrections are taken into account. 

The n — d doublet scattering length 2 a n d is correlated, to some extent, to the 
A = 3 binding energy through the so-called Phillips line [13^ fT4"] . However the 
presence of TNFs could break this correlation. Therefore 2 a n d can be used as an 
independent observable to evaluate the capability of the interaction models to de- 
scribe the low energy region. In Ref. [15] results for different combinations of NN 
interactions plus TNF models are given. We report the results for the quantities 
of interest in Table I. From the table, we can observe that the models are not able 
to describe simultaneously the A = 3, 4 binding energies and 2 a n d- In Ref. [16] 
a comparative study of the aforementioned TNF models has been performed. 
The AV18 [17] was used as the reference NN interaction and the three-nucleon 
interaction models were added to it. Different parametrizations of the URIX, 
TM' and N2LOL TNF have been constructed in order to reproduce, in conjunc- 
tion with the AV18 interaction, £>( 3 H), i3( 4 He) and 2 a n d- In a second step some 
polarization observables in p — d scattering at E[ a b = 3 MeV have been studied. 
In the case of the vector analyzing powers, it was observed that the predictions 
of the different parametrizations appear in narrow bands with different positions 
for each model. Compared to the original AV18+URIX model, the results ob- 
tained using the parametrizations of the N2LOL model were slightly better, in 
particular for A y and %T\\. The results obtained using the parametrizations of 
the TM' were of the same quality. Conversely, the proposed parametrizations 
for the URIX model produced a much worse description of A y and iT%\ than 
the original URIX model. A possible explanation for this fact is the particular 
behavior of the profile functions Y(r) and T{r) used in the construction of the 
model. To this end, in the present paper, we study a different functional form of 
the profile functions Y(r) and T(r). 

2 Three Nucleon Force Models 

In Ref. [15] the description of bound states and zero-energy states for A = 
3,4 has been reviewed in the context of the HH method. In Table Q] we report 
results for the triton and 4 He binding energies as well as for the doublet n — d 
scattering length 2 a n d using the AV18 and the N3LO-Idaho [18] NN potentials 
and using the following combinations of two- and three-nucleon interactions: 
AV18+URIX, AV18+TM' and N3LO-Idaho+N2LOL. The results are compared 
to the experimental values of the binding energies and 2 a n d [19] . 

From the table we observe that the results obtained using an interaction 
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model that includes a TNF are close to the corresponding experimental values. 
In the case of the AV18+TM', the strength of the TM' potential has been fixed 
to reproduce the 4 He binding energy and the triton binding energy is slightly 
underpredicted. Conversely, the strength of the URIX potential has been fixed 
to reproduce the triton binding energy, giving too much binding for 4 He. The 
strength of the N2LOL potential has been fixed to reproduce simultaneously the 
triton and the 4 He binding energies. In the three cases the predictions for the 
doublet scattering length are not in agreement with the experimental value. 

Table 1. The triton and 4 He binding energies B (MeV), and doublet scattering length 2 <w (fm) 
calculated using the AV18 and the N3LO-Idaho two-nucleon potentials, and the AV18+URIX, 
AV18+TM' and N3LO-Idaho+N2LOL two- and three-nucleon interactions. The experimental 
values are given in the last row. 



Potential 


5( 3 H) 


5( 4 He) 


2 

a n d 


AV18 


7.624 


24.22 


1.258 


N3LO-Idaho 


7.854 


25.38 


1.100 


AV18+TM' 


8.440 


28.31 


0.623 


AV18+URIX 


8.479 


28.48 


0.578 


N3LO-Idaho+N2LOL 


8.474 


28.37 


0.675 


Exp. 


8.482 


28.30 


0.645±0.003±0.007 



Following Ref. |16j we give a brief description of the TNF models. Starting 
from the general form 

W= Y, W(i,j,k) , (1) 

i<j<k 

a generic term can be decomposed as 

W(l, 2, 3) = aW a (l, 2, 3)+&W 6 (l, 2, 3)+dW d (l, 2, 3)+c D W D (l, 2, 3)+c E W E (l, 2, 3) . 

(2) 

Each term corresponds to a different mechanism and has a different operatorial 
structure. The specific form of these three terms in configuration space is the 
following: 

W (l,2,3) = Wo(ti •t 2 )(<t 1 -r 31 )(a 2 • r 23 )y(r 3 i)y(r 23 ) 
W 6 (l,2,3) = W Q (n ■ t 2 )[{<t 1 ■ a 2 )y{r 3l )y{r 23 ) 

+ (o-i • r 31 )(cr 2 • r 23 )(r 3 i • r 23 )t(r 31 )t(r 23 ) 

+ (o"i • r 3 i)(<T 2 • r 3 i)t(r 3 i)y(r 23 ) 

+ (<ri • r 23 )(er 2 • r 23 )y(r 31 )t(r 23 )} (3) 
W d (l,2,3) = W (r 3 • n x t 2 )[(<t 3 • <r 2 x ai)y{r 3 i)y{r 23 ) 

+ (<ri • r 31 )(a 2 ■ r 23 )(a 3 ■ r 31 x r 23 )t(r 3 i)t(r 23 ) 
+ (o"l • r 31 )(a 2 ■ r 3 i x cr 3 )t(r 31 )y(r 23 ) 
+ (o~2 ■ r 23 )(a 3 • r 23 x 01 )y(r 3 i)£(r 23 )] , 

with Wq an overall strength. The b- and d-terms are present in the three models 
whereas the a-term is present in the TM' and N2LOL and not in URIX. Here 
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we are interested in the profile functions y(r) and t(r). In the first two models 
these functions are obtained from the following function 

12vr 1 f°° 2 j (qr) 
/0(r) = ^2^./ dqq ^T^/ Aiq) (4) 



where is the pion mass and 



y( r ) = -fo( r ) 



(5) 



t(r) = -y'(r) . 
r 

The cutoff function Fy[ in the TM' or Brazil models is taken as Fy[ = [(A 2 — 
ml) /(A 2 + q 2 )} 2 . In the N2LOL model it is taken as F A = exp(-q 4 /A 4 ). The 
momentum cutoff A is a parameter of the model fixing the scale of the problem 
in momentum space. In the N2LOL, it has been fixed to A = 500 MeV, whereas 
in the TM' model the ratio A/m n has been varied to describe the triton or He 
binding energy at fixed values of the constants a, b and d. In the literature the 
TM' potential has been used many times with typical values around A = 5 m n . 

In the URIX model the radial dependence of the b- and d-terms is given in 
terms of the functions 



Y(r) = e~ x /x^ Y 
T(r) = (l + 3/x + 3/x 2 )Y(r)£ T 



(6) 



with x = m^r and the cutoff functions are defined as £y = £r = (1 — e~ cr ), 
with c = 2.1 fm -2 . This regularization has been used in the AV18 potential as 
well. Since the URIX model has been constructed in conjunction with the AV18 
potential, the use of the same regularization was a choice of consistency. The 
relation between the functions Y{r),T{r) and those of the previous models is: 

Y{r) = y(r) + T{r) 

r 2 (7) 
T (r) = -t(r) . 

With the definition given in Eq.((4]), the asymptotic behavior of the functions 
fo(r), y{r) and t(r) is: 



f (r -> oo) 



3 e 



-X 



m 2 x 



3e~ x ( 1 ' 

j, (r _> oo) _> _ (i + _) (8) 



X* \ X 



3 e~ x / 3 3 

t(r -»• oo) ->• -=■ l + _ + _ 

r z x \ x x z 



To be noticed that with the normalization chosen for /o, the functions Y and T 
defined from y and t and those ones defined in the URIX model coincide at large 
separation distances. 
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The last two terms in Eq.([2]) correspond to a two-nucleon (2N) contact term 
with a pion emitted or absorbed (D-term) and to a three-nucleon (3N) con- 
tact interaction (_E7-term). Their local form, in configuration space, derived in 
Ref. [10J, is 

Wb(l,2,3) = W?(n • r 2 ){( < Ti • a 2 )[y(r 31 )Z (r 23 ) + Z (r 31 )y(r 23 )\ 

+ (<ri • r 3 i)(er 2 • r 3 i)t(r3i)Z (r 2 3) , g . 
+ (fi • r 23 )(a 2 ■ r 23 )Z (r 31 )t(r 23 )} 

WeQ-,2,3) = W(f (n • r 2 )Z (r 31 )Z (r 23 ) . 

The constants W^P and Wq fix the strength of these terms. In the case of the 
URIX model the .D-term is absent whereas the 22-term is present without the 
isospin operatorial structure and it has been included as purely phenomenologi- 
cal, without justifying its form from a particular exchange mechanism. Its radial 
dependence has been taken as Z${r) = T 2 (r). In the N2LOL model, the function 
Zq(t) is defined as 

127T 1 Z 100 

Z (r) = — —5 / dqq 2 j (qr)F A (q) (10) 

with the same cutoff function used before, Fj\{q) = exp(— q A / A ). In the TM' 
model the D- and £^-terms are absent. 

In order to analyze the different short range structure of the TNF models, in 
Fig. [T]we compare the dimensionless functions Zo(r), y(r) and T(r) for the three 
models under consideration. In the TM' model using the definition of Eq. (|10p 
and using the corresponding cutoff function we can define: 



(ii) 

This function is shown in the first panel of Fig. [1] as a dashed line. From the 
figure we can see that, in the case of the URIX model, the functions Zq(t) and 
y(r) go to zero as r — > 0. This is not the case for the other two models and is 
a consequence of the choice to regularize the Y and T functions adopted in the 
URIX. The function Z™(r) has been introduced in Ref. [16] to add a repulsive 
term to the TM' model. In fact, it was shown that without it, the AV18+TM' 
model was unable to reproduce simultaneously the triton binding energy and the 
doublet scattering length for reasonable values of the TM' strength parameters. 



3 Parametrization of the profile functions Zq(t), y{r) and T(r) 



In this section we study a possible parametrization of the profile functions Zq(v), 
y(r) and T(r). The function Zq is defined in Eq. ljlOp . its behavior for small values 
of r can be derived from the expansion of the Bessel function as 



Z (r -> 0) 



127T 1 

~J2^ 



tir 



dqq 



q 2 r 2 



F A (q) 



(12) 



5 5 

r [fm] r [fm] r [fm] 



Figure 1. The Zo(r), y(r) and T(r) functions as functions of the interparticle distance r for 
the URIX (solid line), TM' (dashed line) and N2LOL (dotted line) models. The dotte-dashed 
line shows the one-parameter functions defined in Eqs.f [14)) and ([16}. 



For a sharp cutoff this integral can be approximated as 



12 1 

Z r^O )« — — 
mi 2tt 



A 3 



,A 5 



3 r 30 + " ' 



A 3 



7r \ m 



l A 2 



10 



(13) 



Therefore in the following we propose the one-parameter form of the function 

Z (r) = Z (0)e-^ 2 / 10 (14) 



where Zq(0) can be taken from Eq. fllOp using a particular form of cutoff. In the 

case of sharp cutoff it results: Zq(0) = ^ 

The functions y(r) and t(r) are defined in Eq.([5]). Their short range behavior, 
after expanding the corresponding Bessel function, are 



y(r 0) 



t(r 0) 



12vr 1 1 
m 3 2-7T 2 r 
1 



dqq s 



qr 



(1 



10 



+ 



y(0) + -r 2 t(0) + . . . 



12vr 1 1 
I 2tt 2 ^2 



in 



dqq 



q 2 r 2 



(1 



t(0) 



1 



-r% + ... 



14 



g 2 + m 2 



+ ... 



F A {q) 
q 2 + m 2 



(15) 



where we have introduced the corresponding values at r = and the quantity £2 
in the second term of t(r). In Eq.([8]) the asymptotic behavior of the functions y(r) 
and t(r) are given. Recalling that 3T(r) = r 2 t(r) and considering the short-range 
behavior indicated above, we will analyze the following one-parameter r-space 
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form of the functions y and T 

»M = -^(i + |)(i-.-^) 

(16) 

T(r) = ( i + 3 + 3 \ / _ ^,^(0)1/9^) 

x \ x x l ) V / 

These functions are shown in Fig. [1] with the dot-dashed line. They have been 
calculated using the cutoff of the N2LOL model. As expected they are close to 
the profile functions of the N2LOL potential. Calculations in the A = 3 systems 
using these profile r-space functions are analyzed in the next section. 

4 Results in the A = 3 system 

In the previous section we have presented the one-parameter profile functions, 
Zq(t), y(r) and T(r), obtained from a regularization of their asymptotic form. 
The regularization was perfomed in order to match the short range behavior of 
the profile functions defined in the N2LOL potential. As it is shown explicitly 
in Fig. [TJ these functions are very different from those defined in the Urbana 
potential in which a different regularization was used. In order to study the 
sensitivity in the description of the A = 3 system to different profile functions, 
we construct a modification of the Urbana model in which the profile functions 
defined in Eqs. (|14p and (|16p are used in place of the original ones. The parameters 
^o(O), y(0) and t{0) are calculated using the cutoff of the N2LOL model, Fa = 
exp(— q 4 /A^), with A = 500 MeV. As mentioned, in the URIX model, only the 
&-, d and .E-terms are included. The corresponding strengths are fixed by the 
constants A^^ , and Ar. Their original values are shown in the first row 

of Table (2) Changing the form of the profile functions the values of the constants 
has to be fixed. Three sets of values, selected to reproduce the triton binding 
energy and 2 a n d, are given in Table [2] together with some mean values calculated 
from the triton wave function. 

Table 2. Mean values of the triton kinetic energy, the two-nucleon potential energy V(2N), 
and the attractive, Va(SN), and repulsive, Vr(3N), contributions of the TNF to the triton 
binding energy using the AV18+URIX potential for the specified values of the parameters. In 
the last column 2 a n d is given. 





Ztv 


Ar 


T 


V{2N) 


V A (3N) 


V R (3N) 


Q"nd 


[MeV] 




[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[fm] 


-0.0293 


0.25 


0.0048 


51.259 


-58.606 


-1.126 


1.000 


0.578 


-0.1200 


0.25 


0.0108 


50.110 


-57.360 


-1.747 


0.525 


0.643 


-0.1200 


0.50 


0.0155 


50.193 


-57.328 


-2.097 


0.759 


0.645 


-0.1200 


0.75 


0.0229 


50.331 


-57.211 


-2.735 


1.143 


0.644 



From the table we can observe that with the proposed parametrization the 
attractive part of the TNF is bigger and, due to the fact that the profile functions 
are smoother, there is a reduction of the mean value of the kinetic energy. In 
order to extend further the analysis, p — d scattering observables at Eiab = 3 
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MeV have been calculated using the three sets of constants and the new form 
of the profile functions. The results for the differential cros section, the vector 
analyzing powers A y and iT\\ and the tensor analyzing powers T20, I21 and 
T22, are shown in Figj2] and compared to the predictions of the original URIX 
model and the experimental data. From the figure we can observe that, besides 
a small improvement in A y and iTn, all the models describe the data with 
similar quality. However only the models with the profile functions of Eqs. (114ft 
and ()16p reproduce the experimental value of 2 a n ^. The fact that the vector 
analyzing powers improve very little with the new parametrization can be taken 
as a further evidence that the spin-isospin structure of the URIX is incomplete 
and different forms could in principle be included |21j . 

5 Conclusions 

Due to the fact that some of the widely used TNF models do not reproduce 
simultaneously the triton and 4 He binding energies and the n — d doublet scat- 
tering length, possible modifications of their parametrizations have been ana- 
lyzed. To this end we have used the AV18 as the reference NN interaction and 
we have analyze possible modifications of the URIX model. We have modified 
the regularization of the profile functions Y(r) and T(r) at the origin and we 
have introduced the Zq{t) function in the central repulsive I?-term. We have 
used one-parameter functions that have been chosen to match the short-range 
behavior of the corresponding functions in the N2LOL model. Furthermore the 
strengths of the b-, <i-terms and .E-terms have been fixed to reproduce the tri- 
ton binding energy and 2 aw- Then the predictions for some selected scattering 
observables in p — d scattering at 3 MeV have been compared to the results of 
the original model and the experimental data. We can observe that the descrip- 
tion using the new parametrizations has the same quality of the original model. 
However, with the proposed parametrizations, the AV18+URIX model describes 
correctly S( 3 H) and 2 a n d- This analysis can be consider as a preliminary step in 
a study directed to determine the parametrizations of the profile functions inside 
the three-body force from the experimental data. Investigations in this direction 
are underway. 
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